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Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map 
the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity (Bo) and curve 
and the distance to the Galactic center (Ro). Here, we present an update on our first results, implementing a recent change 
in the knowledge about the Solar motion. It seems unavoidable that the lAU recommended values for Ro and Go need 
a substantial revision. In particular the combination of 8.5 kpc and 220 km s~^ can be ruled out with high confidence. 
Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr A* 
gives best values of Ro = 8.3 ± 0.23 kpc and Go = 239 or 246 ± 7 km s~\ for Solar motions of V© = 12.23 and 5.25 
km s~^, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) 
Survey. 
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1 Introduction 



The Milky Way is a barred spiral galaxy, as seen from ob- 



servations of CO and Hi gas (e.g. B urton 198 8i|Dame et al. 
|2001| ), and star counts (e.g. |Benjaniin et al.||2005| ). How^ 
ever, our location in the Galaxy makes it difficult to deter- 
mine the number and positions of spiral arms, the length 
and position of the central bar, and the rotation curve. As a 
result, even the most fundamental parameters of the Milky 
Way, such as the distance to the Galactic center Rq and the 
rotation speed 6o are still not known with high accuracy. 
However, these values are not only important for Galac- 
tic astronomy, but also for a wide range of different fields. 
This includes the interpretation of the proper motions of 



the Magellanic Clouds (Besla et al. |2010[ 


Diaz & Bekki 


2011 


|Peebles||2010l|Ruzicka et al.||2010t|S 


lattow & Loeb 


2009 


) and galaxies in the Andromeda subgroup ( Brunthaler 


etal. 


20051 120071 Ivan der Marel & Guhathakurta 2008 ), the 



motion of the Sun relative to the cosmic microwave back- 



Corresponding author: e-mail: brunthal@mpifr-bonn.mpg.de 



ground (e.g.| Loeb & Narayan|20Q8 ), and even on the inter 
pretation of dark matter direct detection experiments ([Foot 



,2010,McCabe 2010 ). 

In recent years, many large scale surveys have covered 
the Galactic Plane in all wave bands from radio to gamma 
rays. However, almost all of these surveys are all only two 
dimensional, and using them to construct a three dimen- 
sional model of the Milky Way is not trivial, mainly due to 



large uncertainties in distance measurements ( G eorgelin & 
Georgelin|1976||Hou et al.|20Q9| ). These also affect the in 



terpretation of these surveys since most astrophysical quan- 
tities, such as linear size, mass, and luminosities, strongly 
depend on the distance to the object. 

2 Galactic distances 

In Galactic astronomy, the kinematic distance is very com- 
monly used. Here, a distance is deduced from a measured 
line-of- sight velocity, a rotation model of Milky Way, and 
the assumption that the object has no peculiar motion. Apart 
from a distance ambiguity in the inner Galaxy (near and far 



© 2010 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



790 



Brunthaler, Reid, Menten et al.: Mapping the Milky Way 



kinematic distance), errors in the rotation model and pecu- 
liar motions can lead to very large errors in the kinematic 
distance. One extreme example is the star forming region 
G9.62+0.20. It is located at a distance 5.2 ± 0.6 kpc ( |Sanna| 
[et al. 2009 ) , while the near and far kinematic distances place 
the source at 0.5 and 16 kpc, respectively. This source shows 
a very large peculiar motion which may be induced by the 
central bar and is responsible for the large error in the kine- 
matic distance. 

Certainly, the trigonometric parallax is the most funda- 
mental method, since it is based on pure geometry without 
any astrophysical assumptions except the well known orbit 
of the Earth about the Sun. However, the measurement of 
an accurate trigonometric parallax requires extremely high 
astrometric precision. Friedrich Wilhelm Bessel was able to 
measure the first stellar parallax of the star 61 Cygni ( |Bessel| 
1838a|b ). A huge step forward came with the ESA's Hip- 



Source 



Parallax 
[//as] 



Reference 



parcos satellite (Ferryman et al. 1997). It provided astro- 
metric accuracies of the order of 1 milliarcsecond, which 
allows distance estimates in the Solar neighborhood out to 
100 pc with 10% accuracy, i.e. a very small portion of the 
Milky Way. Following the footsteps of Hipparcos, ESA's 
new Gaia mission ( jPerryman et al.||2001 ) will be launched 
in late 2012 to provide astrometry of up to 1 billion stars 
in the Milky Way with parallax accuracies up to 7 /ias if 
it achieves specifications, a factor 100 better than Hippar- 
cos. Although Gaia will revolutionize our view of the Milky 
Way, the strong extinction by gas and dust in the Galac- 
tic plane and in particular the spiral arms will prevent Gaia 
from making significant progress on the spiral structure of 
the Milky Way. Since radio waves are not affected by inter- 
stellar extinction, radio astronomy can come to the rescue 
and fill the gaps. 

3 VLBI Astrometry 

Recently developed calibration techniques for Very Long 
Baseline Interferometry (VLBI) have improved the accu- 
racy of astrometric VLBI observations significantly. When 
these techniques are applied to VLBI networks like the 
NRAO Very Long Baseline Array (VLB A) in the US, the 
European VLBI Network (EVN) in Europe and China, or 
the VLBI Exploration of Radio Astrometry (VERA) array 
in Japan, it is now possible to measure parallaxes of radio 
sources with accuracies of the order of 10 /ias, comparable 
to the expected accuracy of Gaia. 

Possible target objects for VLBI astrometry are either 
radio continuum sources or strong maser emission. Most 
radio stars are relatively weak, and the sensitivity of cur- 
rent instruments limits the detection to sources within a few 



hundred parsec (e.g.|Loinard et al. 2007[ Menten et al. 2007 



Torres et al .[2007[). However, sensitivity upgrades of exist- 
ing VLBI instruments (e.g. [Ulvestad et al.|2010 ) and even- 
tually the Square Kilometer Array (*Fo malont & Reid|2004| ) 
will extend the accessible range to even weaker and more 
distant sources. 



W3(0H) 512 ± 10 

Orion Nebula 2425 ± 35 

S269 189 ±16 

VYCMa 876 ±76 

NGC281 355 ±30 
S 252 476 ± 6 

G232.6+L0 596 ± 35 

CepA 1430 ±80 

NGC7538 378 ±17 

W51IRS2 195 ±71 

G59.7+0.1 463 ±20 

G35.2-0.7 456 ± 45 

G35.2-1.7 306 ±45 

G23.0-0.4 218 ± 17 

G23.4-0.2 170 ± 32 

G23.6-0.1 313 ± 39 
WB 89-437 167 ± 6 

IRAS00420+5530 470 ± 20 

Table 1 List of 1 8 sources which were used to fit a new 



Xu et al. ( 2006 ) 
Menten et al. ( 2007j) 
Honma et al. (20^ 
Choi et al. ( 2008| 
Sato et al. (2008 ) 
Reid et al. ( 2009a I 

Reid et al. ( 2009a I 

Moscadelli et al. ( 2009 
Moscadelli et al. (200? 
Xu et al. ( 2009 ) 



Xuetal. (2009 ) 
ZEanget al. (200|) 

Zhang et al. (2005| 

iBrunthaler et al. ( 2009 1 
Smnth aler e t al..(2009 1 
Bartkiewicz 



|Hachisuka 
Moellenbrock et al. (2009 ( 



model of the Galaxy in |Reid et"aL]p009b] ). 



Molecules with maser lines suitable for astrometry are 
hydroxyl (OH), methanol (CH3OH), water (H2O), and sil- 
icon monoxide (SiO). Since OH masers at a low frequency 
of 1.6 GHz are strongly affected by interstellar scattering 
and ionospheric effects, they are not the optimal targets for 
VLBI astrometry. SiO masers at 43 GHz are less affected 
by scattering and the ionosphere, but these observations are 
more vulnerable to tropospheric phase errors and require ex- 
cellent weather conditions. Nevertheless, both masers have 
been used to measure accurate distances to evolved stars 
(e.g. |Choi et al.|2008[[van Langevelde et al.pOOO] ). 

The most suitable targets are methanol (6.7 and 12.2 
GHz) and water (22 GHz) masers. They are very strong 
and are found in high mass star forming regions (HMS- 
FRs) mainly in the spiral arms of the Galaxy. The catalog 
of 6.7 GHz methanol masers by |Pestalozzi et al.| ([2005 ) lists 
more than 500 sources, and many more were found recently 
with Arecibo ( Pandian et al. 2007), Effelsberg ( |Xu et al.] 
2008 ), and in particular in the methanol multi beam survey 



on the Parkes telescope ( ^Caswell et al.||2010[ |Green et aL 
2010| ). Water masers can be even stronger andJValdettaro 
et al.|pOOT]) list more than 1000 sources in the Galaxy. 



Parallaxes using the 6.7 GHz methanol transition have 
been reported so far only with the EVN ( |Rygl et "aL 
2010a b). while the VLBA can currently observe only the 
12.2 GHz methanol and 22 GHz water maser lines. One ex- 
ample is the distance to W3(0H), a high mass star forming 
region in the Perseus spiral arm of the Milky Way. Two in- 
dependent parallax measurements with the VLBA of bright 
methanol ( |Xu et al.|2006| ) and water ( [Hachisuka et al.|2006| ) 
masers yielded consistent distance estimates of 1.95 ± 0.04 
and 2.04 ± 0.07 kpc. 

Another example that demonstrates the high quality of 
VLBI parallax measurements is the Orion Nebula. Four 
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independent measurements with two different instruments 



give, within their joint errors, consistent results. Sandstrom 



et al.| ( |2007| ) observed one radio star with the VLBA and 
obtained a distance of 389^21 P^' while 



Hirota et al. 



(2007) 



used VERA to measure the parallax of water masers in the 
nebula (437 ±19 pc). The most accurate measurements of 
414 ± 7 pc (Menten et al. 2007, VLBA, 4 radio stars) and 
418 ± 6 pc ( Ki m et al. 2008 , VERA, SiO maser) also agree 
with each other and with the modeling of the orbit of the 



binary ©^OriC from near-infrared interferometry (Kraus 
|etal.|2007| 434±12pc). 



4 A new model for the Milky Way 



4.1 Measurements so far 



In |Reid et al.| ( |2009b] ), we used 18 sources, with published 
parallaxes at that time (see Table [T] and Fig.[T]), to analyze 
the spiral structure and the rotation of the Milky Way. We 
were already able to locate the Local arm and the Perseus 
spiral arm from Galactic longitudes / = 122° - 190°. The 
measured pitch angle of the Perseus spiral arm of 16° ± 
3° favors four rather than two spiral arms for the Galaxy. 
Furthermore, we find that most of the HMSFRs are closer 
than their kinematic distance. Note that the same data set 
has been also analyzed by Bovy et al.|( 2009|), 
Binney| ( |2QT0l ), and |Bobylev & Bajkova[pQlQ 



McMillan & 




Fig. 1 Artist conception of the Milky Way (R. Hurt: 
NASA/JPL-Caltech/SSC) with the positions of the 18 
sources from Table [T] Sources in blue are 12.2 GHz 
methanol masers measured with the VLBA, while the green 
circles indicate H2O masers observed with the VLBA or 
VERA and the Orion Nebula. The red symbols mark the 
position of the Galactic center and the Sun. Error bars are 
shown for all sources, but smaller than the symbol size for 
most sources (taken from ,Reid et al.|2009b| ). 



4.2 Solar Motion 

Astrometric observations not only yield the distance but also 
the proper motion of a source. Together with the known po- 
sition and line-of- sight velocity, this gives the full 6 dimen- 
sional phase space information for the observed sources. 
However, all measurements are heliocentric. Therefore, the 
peculiar motion of the Sun relative to the LSR is needed to 
convert the measured heliocentric into Galactocentric mo- 
tions. For over a decade, the values of \Jq= 10.00 ± 0.36, 
V0= 5.25 ± 0.62, and Wq= 1 Al ± 0.38 km s"^ for the 
peculiar motion of the Sun derived from Hipparcos data by 



Dehnen & Binney (19981 have been widely used. When us- 



ing this Solar motion, we find that the HMSFR s rotate on 
average 15 km s~^ slower than the Milky Way (ReidetaL 
|2009b| ). 

This controversial result has sparked new interest in the 
Solar motion and Binney ( 2010| ) argues that this slower rota- 
tion is partly induced by an erroneous value of VQ.|Mignard 



( |2000| ) and Piskunov et al. ( 2006|) already favored a signifi 
cantly higher value of V© , in the range of ~ 12 km s~ ^ . This 
higher value, which is similar to the pre-Hipparcos value 
(e.g. Mihalas & Binney [1981 1), has been also found in a 
number of recent studies (Co§kunoglu et al.|201 if [Francis &| 
Anderson|20Q9l|Schonrich et al.,2010) , while the values of 
U0 and W0 have not changed considerably (see also|Fuchs| 
et al.|2009 ). Although the issue of the V© component of the 



Solar motion is probably still not solved, it seems likely that 
the value of V© is probably closer to 12 km s~^ than to 5 
km s~^. 

4.3 Revised values 



The analysis in |Reid et aL] ( |2009b| ) was based on the old 
Hipparcos value of the Solar motion from Dehnen & Bin- 
'ney'(fT998]). Here we give an update using the new values of 
Schonrichetal.|(|20T0]), which are U© = 11.10 ± 1.2, Nq= 



12.24 ±2.1, and W0= 7.25 ± 0.6 km s-^ Using this value, 
the HMSFRs rotate 8 ± 2 km s~^ slower than the Milky 
Way. By fitting our measurements to a model of the Galaxy, 
we can also estimate the distance to the Galactic center Rq 
and the circular rotation speed 60. Assuming a flat rotation 
curve, we get values of Rq = 8.4 ± 0.6 kpc and 60 = 247 
±16 km s~^. A similar analysis allowing for different ro- 
tation curves yields results in the range of Rq = 7.9 - 9 kpc 
and 60 = 223 - 280 km s~^ . However, the ratio of these two 
values is much better constrained with a value of 9o/Ro = 
29.4 ± 0.9 km s~^ kpc~^. Clearly, the number of sources 
is currently not large enough to constrain different rotation 
curves. 

The slower average rotation of the HMSFRs together 
with the higher overall rotation of the Milky Way also ex- 
plains why most of the kinematic distances are larger than 
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the true distances. A description that takes into account 
these two effects and gives (in general) more accurate kine- 
matic distances, and a more realistic distance uncertainty 
estimate is also presented in^ Reid et al.| (|2QQ9b). It should 
be noted, however, that these revised kinematic distances 
can be still unreliable, because sources may have very large 
peculiar motions. 

4.4 Independent measurements of Rq and Oq/Rq 

The distance to the Galactic center Rq has been the target 



of numerous investigations over the last decades (see Reid 



1993 for a review). Most measurements in the range be- 
tween 7.5 and 8.5 kpc and direct geometric distance esti- 
mates are rare. The accurate determination of stellar orbits 
in the Galactic center yields now consistent values from 
two groups of 8.4 ± 0.4 kp c (|Ghez et al.|[2QQ8 ) and 8.33 
±0.35 kpc ( [Gillessen et al.||2009| ). A trigonometric paral- 
lax to water masers in SgrB2, a star forming region located 
within 150 pc from the Galactic center, also yields a consis- 
tent value of Rq = 7.9t.o\Y ^P^' ^i^^ a larger uncertainty 
( |Reid et al.|2009 c'). Further observations of the water masers 
in SgrB2 will improve the accuracy to a value comparable 
or even better than from the stellar orbits. 

The ratio of 6o and Rq is known with better than 



1% accuracy from the proper motion of Sgr A* ( |Reid &| 
Brunthaler 2004| ). This motion in the Galactic plane of 
6.379 ± 0.026 mas yr~^ is a combination of Galactic ro- 
tation and Solar motion, and corresponds to a value for 
(9o+Vq)/ Ro of 30.24 km s~^ kpc ~^ Using the Solar mo- 
tion from I Schonrich ( 2010 1), and assuming the dis- 
tance of 8.4 kpc, yields a value of Oq/Ro = 28.79 ± 0.26 
km s~^ kpc~^, where the uncertainty is dominated by the 
uncertainty in the Solar motion. This value depends only 
weakly on the assumed Rq (in the correction of the So- 
lar motion), and does not change by more than 0.7%, even 
for extreme values of Rq of 7.5 or 9.5 kpc. This value is 
perfectly consistent with our results from the independent 
maser parallaxes but significantly larger (>11 a) than the 
lAU value of 25.88 from the combination of Rq = 8.5 kpc 
and Oo = 220 km s~^. Hence, one requires an Rq of less 
than 7.7 kpc to have a rotation velocity of 220 km s~^. 

4.5 Conclusion and best values for Rq and Oq 

The measured maser parallaxes, and the combination of the 
stellar orbits with the proper motion of Sgr A* provide inde- 
pendent and consistent evidence for a higher rotation veloc- 
ity of the Milky Way. Therefore, it seems unavoidable that 
the lAU recommended values for Rq and 6o need a sub- 
stantial revision. In particular the combination of 8.5 kpc 
and 220 km s~^ can be ruled out with high confidence. 

The weighted average of the four direct measurements 
for Rn of 8.4 zb 0.4 kp c fGhez et al. 20 08]), 8.33 ± 0.35 kpc 



( [Gillessen et al.|2009|, 7 . 910 7 ^^ P^^ ( ^Reid et al.|2009c| ), and 
8.4 zb 0.6 kpc |Reid et al.] ( |2009b| ) gives best value of 
Ro = 8.3 ±0.23 kpc. 



This translates then into a rotation speed of 

Go = Oo/Ro X Ro = 28.79 km s"^ kpc"^ x 8.3 kpc 
= 239 ±7 km s"^ 

These values are almost identical to values presented by 
McMillan pOll I and marginally consistent with the rota 
tion speed of 224 ±13 km s~^ provided by Koposov et al. 



2010| ) from modelling the stellar stream GD-1. For the old 



Solar motion of V© = 5.25 km s~^, 6o would increase to 
246 ± 7 km s~^. The uncertainties are formal uncertainties 
and are dominated by the uncertainty in Ro . The true uncer- 
tainties are probably a factor of V2 larger, since the two stel- 
lar orbit distances may have similar systematic errors (they 
are independent measurements with different instruments, 
but of the same source). 

A higher value of 9o is also supported by the observa- 
tion of a retrograde-rotating and metal-poor component in 
the stellar halo of the Milky Way by peason et al.| ( |2010| ) 
when assuming a value of 6o = 220 km s~^. A value of 
6o ~ 240 km s~^ would explain this apparent retrograde 
rotation. 

5 The Bar and Spiral Structure Legacy 
Survey 



Bess 



Bar and Spiral Structure Legacy Survey, 
a VLBA Key Science Project 



Motivated by these very encouraging results from only 
18 sources, we have started a much larger project, the Bar 
and Spiral Structure Legacy (BeSSeL|^ Survey, a VLBA 
Key science project. The goal of BeSSeL, named in honor 
of Friedrich Willhelm Bessel who measured the first stellar 
parallax, is to measure accurate distances and proper mo- 
tions of up to 400 high mass star forming regions in the 
Milky Way between 2010 and 2015. This will result in a 
catalog of accurate distances to most Galactic high mass 
star forming regions visible from the northern hemisphere 
and very accurate measurements of fundamental parameters 
such as the distance to the Galactic center (Rq), the rotation 
velocity of the Milky Way (6o), and the rotation curve of 
the Milky Way. Additionally, maps of the maser distribu- 
tion from the first epoch of each source will be published on 
the BeSSeL website shortly after the observations. 



http://www.mpifr-bonn.mpg.de/staff/abrunthaler/BeSSeL/index.shtml 
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Fig. 2 Similar to Fig. [jj but showing all sources currently measured (green), including unpublished sources, and all 
sources observed in the first year of BeSSeL (red), based on their kinematic distances. 



The BeSSeL Survey will first target 12.2 GHz methanol 
and 22 GHz water masers. Once the VLB A is equipped 
with new receivers that also cover the 6.7 GHz methanol 
maser line (presumably in 2012) these masers will be also 
observed. In early 2010, preparatory surveys started with the 
Very Large Array to obtain accurate positions of the target 
water masers, and with the VLB A to search for extragalactic 
background sources near the target maser sources (Immer 
et al.|201l] ). The first parallax observations started in March 
2010, and first results are expected in mid 201 1 (see Fig.|2]). 
In parallel, the VERA array will observe additional H2O 
and SiO masers throughout the Galaxy. Combined with 



complementary efforts in the southern hemisphere with the 
Australian Long Baseline Array, this will result into a de- 
tailed and accurate map of the spiral structure of the Milky 
Way. The superior sensitivity and the large field-of-view of 
the Square Kilometer Array, which will also cover the 6.7 
GHz methanol maser line, will reach even higher astromet- 
ric accuracies for methanol masers, radio continuum stars, 
and pulsars. This will result in a very detailed map of the 
the spiral stucture in the southern hemisphere. 
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